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The possible dependence of fundamental couplings and mass ratios on the gravitational potential has 
been bounded by comparing atomic clock frequencies over Earth's elliptical orbit. Here we evaluate 
bounds on such dependence from Eotvos-type experiments that test the Weak Equivalence Principle, 
including previously neglected contributions from nuclear binding energy. We find that variations 
of fundamental parameters correlated with the gravitational potential are limited at 10~*-10~^, an 
improvement of 2-3 orders of magnitude over atomic clock bounds. 
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Both the Standard Model of particle physics and the 
theory of General Relativity are constructed on the as- 
sumption of Local Position Invariance (LPI): that locally 
measurable properties of matter do not vary over space 
and time. LPI forms part of Einstein's equivalence prin- 
ciple. While most tests of LPI concerning particle cou- 
pling strengths and masses have yielded results consistent 
with zero variation there is a significant indication 
of cosmological variation in the fine structure constant 
a deduced from astrophysical absorption spectra 0, [1] . 
Such claims are controversial [1, H, [y, 0] and await im- 
provements in data and analysis. The proton-electron 
mass ratio /i = mp/me has also been probed: analysis of 
molecular Hospectra indicates a moderately significant 
variation [i,l9[ at high redshift, while recently a strong 
bound on variation at lower redshift has been derived 
from the NII3 inversion spectrum [lo| . 

It is also of interest to probe possible variations at 
the present time within the Solar System. Indeed any 
theory with an underlying Lorentz invariance, in which 
time variation may occur, should also allow spatial varia- 
tion. Recently bounds have been set 11 , 13] on variations 



of fundamental parameters correlated with the gravita- 
tional potential U ^ by comparing atomic clock frequen- 
cies over several months while the Earth moves through 
the Sun's gravitational field 1^ For independently 



varying parameters Gi we may define couplings ki via 



AlnG, = fc,AC/ 



(1) 



for small changes in U ; the annual variation due to 
Earth's orbital eccentricity is AJ7 ~ 3 x 10~^°. The 
couplings of a, /i and the light quark mass mg/Ac = 
{ruu + md)/2Ac, where to„ and nid are the up and down 
quark current masses and is the invariant strong in- 
teraction scale of QCD [l^ , were found to be consistent 
with zero [TTl, with uncertainties 



{a(fc„), a(fc^), aikg)} = {3.1, 17, 27} x 10" 



(2) 



Now as pointed out in [16|, a spatial gradient of cou- 
plings or masses will lead to anomalous accelerations of 
test bodies, due to the dependence of their mass-energy 
on the varying parameters. With the gravitational ac- 
celeration g given by — VJ7, we have VlnG^ — ^kig. 



Rather than moving on a geodesic with acceleration g, 
a freely-falling body of mass M{Gi{x)) will experience 
additional acceleration a = — ( VM) /M [l3| , as if moving 
in a potential V{x) — M{Gi{x)). 

Units and definitions For a test body with a 
given space-time dependent mass Mh(x) = A/o(l + 
AAf(x)/Afo), we may perform a redefinition of units 
[lit such that the value of A/;, is everywhere constant. 
The Newton "constant" Gn is then in general space-time 
dependent, and the body moves along geodesies of the 
rescaled metric = (Mo(x)/A/o)^5^i, [ig. In the New- 
tonian Hmit we have 500 - -1 - 2U{x) - 2AM{x)/Mo, 
and the previous result for a is recovered. For two test 
bodies of different x-dependence such a redefinition can- 
not be performed and the differential acceleration db — dc, 
violating the Weak Equivalence Principle (WEP) or uni- 
versality of free fall, is a physical signal probing spatial 
gradients of Gi. 

For convenience we will express masses in Planck units 
where Gn is constant. The acceleration is then Il9[ 



EdlnMb^. ^ainAffc 
— — -VlnG. = )^^,^A:,5, 



ainG,- 



(3) 



and the differential acceleration is 



Vb- 



\ab 



\9\ 



E 



ain(Afb/Af,), 



91nG, 



(4) 

which defines the sensitivity coefficients A J Each ex- 
perimental limit on rj then bounds some linear combina- 
tion of the couplings ki. To bound n independent cou- 
plings ki , at least n independent pairs of test bodies with 
different derivatives d In M^''^/ d In Gi are needed. 

We may also change basis from one set of parameters 
Gi to another GL such that ^91 



Fikki 



F,i 



ainG. 
dlnGi 



(5) 



for couplings fc^ in the new basis. 

It is also essential to consider a complete set of param- 
eters Gi which account for all non-negligible sources of 
variation in Mb/Mc, otherwise measurements of rj may 



2 



be incorrectly interpreted. Couplings to baryon number 
A, lepton number L = Z and electromagnetic self-energy 
Ef.y ,, ~ acZ{Z—\)IA^^'^ have commonly been considered 
[2C| . However, there are other contributions to A/^ which 
depend on different, linearly independent functions of A 
and Z. These arise from the strong nuclear binding en- 
ergy, which for large A is well approximated by the liquid 
drop model: 



Br, 



A 



ay 



as 



-aA- 



{A~2Zf 



-ap- 



(6) 

where the coefficients of volume, surface, asymmetry, 
and pairing energy respectively are {av,as,aA,ap} ~ 
{15.7, 17.8, 23.7, 11.2} MeV 21]. Such terms could a pri- 
ori couple to U and induce differential acceleration. 

Statistical treatment We may construct [2Q . i22j | a like- 
lihood function 



{b-c} 



<Jb-c 



(7) 



where the sum over {b — c} counts independent exper- 
imental results with central values rjh-c and la uncer- 
tainties (Jfj-c- The likelihood C is proportional to e^^ 1"^ . 
A "metric" gij may be defined by the coefficient of the 
quadratic term kikj in Eq. ([7|: 



{b-c) 



— 2 A b—c \ b—c 
b-c-^i 



The best-fit point with minimal — Xmin then 

h.min = (Q^^)t, V crr_l?76-cA-"'= (9) 



{b-c} 



and the excess of above this minimum is given by 

Ax' = - xLn = 9^J^h^k,, (10) 

where Aki = ki — ki^min- Error ellipsoids around the cen- 
tral values are then surfaces of constant Ax' . The eigen- 
vectors of g give the linear combinations of ki pointing 
along the principal axes of the ellipsoids. 

We may project the likelihood onto each ki separately 



111 |23| | by integrating over kj {j ^ i), resulting in a 



Gaussian distribution with variance 

2 Ma 
cr. 



det g ' 



where Ma is the minor of the element ga. On changing 
the basis of parameters via Eq. ([5]) , the metric becomes 



fffcr, 



Nuclear parameters and couplings To calculate A^"*^ 
we identify the parameters controlling possible variations 
in mass ratios. We first consider "nuclear parameters" 
which characterize physics at low energy. 

The dimensionless parameters Xj and their couplings 
to U are: the fine structure constant a (coupling ka)] the 
electron mass me/rnjsi (fce); the nucleon mass difference 
Sn/ttin (ksN)', and the nuclear surface tension as/rriN 
{kas)i where Sn = nin—mp and m^r = (m„-fmp)/2. Fur- 
ther couplings may be defined for other terms in Eq. ([6|). 
The normalized mass per nucleon of a body is then 



AruN \ 2 / niN riiN 



Z{Z-l) ac 



ay 



.4/3 • — 

A"*'-^ TUN riiN rriN 



(13) 



where fp = Z/A, the coefficient ac — 0.71 MeV varies 
proportional to a, and the remaining terms are sublead- 
ing in A. Expanding in small quantities we find 

Aln^ = Aln^ 

-i^Alni^ + ^Aln^) A,^J, 



A In a Vc^^|n^ + ^ A In ^ A.-.A" ^3, 



(14) 



(8) where A;,_c denotes the difference of a quantity between 
two test body materials The sensitivity coefficients 
Aj~'^ may then be read off. We cannot distinguish be- 
tween couplings to Sn/itin and rUe/rnN, since their coef- 
ficients are both proportional to At-cfp- Thus we define 
Qn = S]\[ — me (the kinetic energy in neutron decay) with 
a coupling kgn = A\n{Qn/miq) / AU . Its sensitivity co- 
efficient is then \^q^ = -{Qn/mN)Ab-cfp- 

Results We summarize current Eotvos experiments, 
specifying their pairs of test materials. 

Schlamminger et al. [il: Be-Ti, 77 = (0.3 ± 1.8) x 
10^^3 Be-Al data were also taken (to be published). 

BaeBler et al. '2^: Fe-Si02, rj = (0.5 ± 9.4) x 
10^ ^3 Test bodies were 'Earth core' Fe-Cr alloy and 
'Moon/mantle' silica body with a small fraction of Mg. 
Y. Su et al. [2^: 

Be-Al, 7] = (-0.2 ± 2.8) x IQ-^^ 
Be-Cu, T] = (-1.9 ± 2.5) x IQ-^' 
Si/Al-Cu, 7? = (5.1 ± 6.7) X 10"^' 

(11) (Si dominant in Si/Al body). 

Braginsky and Panov [27*1 : Pt~Al, 77 = (-0.3±0.4)x 
10^ Here we rederived a 68% confidence interval from 
the 7 quoted data points. The meaning of the sign is 
ambiguous; we treat the result conservatively as a limit 

(12) < 7 X 10-13 (la). 
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The relevant functions of A and Z are listed in Table |T] 
(4l| . The resulting sensitivity metric g in the basis fc/ = 



^0.00015 -0.0020 0.0101 ' 
-0.0020 0.8696 -2.030 
y 0.0101 -2.030 4.949 



X 10 



20 



(15) 



The eigenvalues are = (7.0 x lO^-*, 3.2 x 10^^, 5.8 x 
10^"), and we find that the linear combinations 



/0.998 
0.064 0.923 
y0.002 -0.382 



0.058 -0.026^ 
0.381 
0.924 



fco 



(16) 



of couplings are bounded at the 10~^, 10^^ and 10^^° 
level respectively "i^ . 

Projecting the likelihood onto each direction fc/ via 
Eq. ()11|) . we obtain null bounds on three couplings with 
uncertainties 

{a(fcQ„), a(fc„), a{kas)} = {38, 2.3, 1.0} x 10"^ (17) 

at least two orders of magnitude stronger than current 
atomic clock bounds, although we note that WEP ex- 
periments probe slightly different linear combinations of 
parameters from clock comparisons ;28j . 

Including asymmetry energy In the above calculation 
we neglected the asymmetry term in Eq. ^ whose sensi- 
tivity coefficient is X^'/ = ^Ab_c((A-2Z)2/A2). If we 
include this contribution with an independent coupling 
kaA, we find a 4-by-4 sensitivity metric g with eigenval- 
ues (3.4 X 10^^, 3.1 X 10l^ 5.6 x 10^®, 5.8 x 10^"), and the 
bounds projected onto each coupling separately are 

{a(fcQ„), a(fc„), (T(fcas), (T(fcaA)} = {52, 14, 5.9, 19}x 10-^. 

(18) 

The increase in uncertainty occurs because four experi- 
ments are now required to solve for the couplings fc/ , but 
there are currently only three independent bounds with 
precision exceeding 10~^^. 



TABLE I: Sensitivity functions for test bodies. For Si02 we 
give a mass- weighted average. 



Material 



Be 
Al 
Si 



SiOa 



Ti 
Fe 
Cu 
Pt 



A 



fp - 0.5 



(A-2ZY 



9 

27 
28.1 



4 
13 
14 



-0.0556 
-0.0185 




0.64 
1.93 
2.14 



0.481 0.0123 
0.333 0.0014 
0.329 0.0000 



21.6 10.8 







1.74 



0.365 0.0000 



47.9 
56 
63.6 



22 
26 
29 



195.1 78 



-0.042 
-0.0357 
-0.044 
-0.1 



2.65 
3.03 
3.20 
5.31 



0.275 0.0066 

0.261 0.0051 

0.251 0.0078 

0.172 0.0402 



Given the common origin of the surface and asymmetry 
terms in nuclear forces, we may assume that the coeffi- 
cients as/mjv and aA/fniq vary in the same way. Then 
we have kaA = kas = fcnuc and we may rewrite with 



^as'^^aS + >\AkaA 



^nuc'^rLuc 



(19) 



where Xnuc = XaS + XaA- We obtain a 3-by-3 metric 
in the basis [kgn, ka, knuc) with small corrections com- 
pared to Eq. The eigenvalues are (6.5 x 10^'*, 5.6 x 
10^*, 6.0 X 10^") and the bounds projected onto each cou- 
pling separately are now 

{crikQn), c7{ko,),a{knuc)} = {39, 1.7, 0.8} x 10~^ (20) 

Hence with reasonable physical assumptions, bounds on 
kj are stable, or even improve, under increasingly precise 
approximations to nuclear binding energy. 

Fundamental parameters and couplings We 
now consider how the "nuclear parameters" 
— {Qn/rnj^,a,as/rnj^) depend on parameters 
Gk of the Standard Model of particle physics. These 
parameters are: the fine structure constant a (cou- 
pling fc^); the electron mass rUe/Ac (fcg); the light 
quark mass m^/Ac {k'^)] the up-down mass difference 
5q/Ac = [rud — m„)/Ac (k'g^); and the strange quark 
mass ms/Ac, where we normalise to the QCD strong 
interaction scale. 

With three independent "nuclear" couplings 
{kQn,ka,knuc), wc Can only bound the couplings of 
three independent combinations of fundamental pa- 
rameters. The dependence of nucleon masses and 
nuclear forces on m^/Ac is subject to large uncertainties 
(possibly 100%): a reliable calculation of the effects of 
varying m^/Ac is a formidable challenge. To proceed 
further, we must assume that m^/Ac does not couple 
significantly to U. 

The dep endence of nucleon masses on Gk was discussed 
in l2 9| . l30l| using results from chiral perturbation theory 
[3l|, |32]. We find, suppressing the term in TOs/Ac, 



Aln 



A, 



0.048 Aln ^- 

A, 



(2.7x10"") Alna- 



(21) 



For the isospin-violating mass difference Qn we derive 



Aln 



-^~2.6Aln^-0.65Aln^ 

m^r Ac Ac 



-0.97Alna-0.048Aln— 2-, (22) 

thus fcQ„ ~ 2.6A:^^ - 0.97fc^ - 0.048fc^, where k'g^ = k'g^ - 
0.25fcc. Clearly also ka — k'^. 

Nuclear binding energy The dependence of strong nu- 
clear binding energy on QCD parameters is subject to 
significant uncertainties. For the deuteron a detailed ef- 
fective field theory was applied H, 11] giving a depen- 
dence on the pion mass m-^ of 

A ln(BD/Ac) = rA ln(m^/Ac) ~ ^ A ln(m,/Ac), (23) 
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where — 10 < r < —6 and we use the leading order depen- 
dence oc ^rriq. Taking this relation as a first approx- 
imation to the dependence of other nuclear binding ener- 
gies [3^ we estimated dBi/dm^ = fi{Ai — l)r(_BD/m-ir), 
with fi 1. This is consistent with estimates from re- 
alistic nuclear interaction models in small nuclei [ssj if 
r ~ —7. An effective theory treatment of the central 
force in larger nuclei yields a strong negative dependence 
of a A on m^r/Ac [36[ which would only strengthen the 
resulting bounds. Thus 



A, 



iZ^ai^li^Aln^ 
2B, A. 



-0.9Aln^, 



'(24) 

taking Bi/{Ai — 1) ~ 8.5 MeV. If, as argued above, all 
coefficients in Eq. ([6]) react in a similar way to variation 
of niql Kc, knuc is determined simply by the behaviour of 
Bi/mpf and we have fc„„c — — 0.9fc^. Thus the matrix 
Fjk giving the dependence of Xj on Gk is 



(25) 



where the columns correspond to variations of parame- 
ters {Sf/Ac, a, rUg/Ac). 

Transforming the sensitivity metric g (including the 
asymmetry energy contribution) via Eq. (|12p we obtain 
a metric g' for the couphngs {k'gj,k'^,k'g). This yields null 
bounds with uncertainties 




{aik'sj),cTik'J,a{k'q)} - {14, 1.7, 0.9} x 10" 



(26) 



This result is subject to theoretical uncertainty in the 
dependence of nuclear properties on quark masses, and 
is obtained with the assumption that variation of rris/ Ac 
is negligible [i^l- However, it is clear that gravitational 
experiments testing WEP provide by far the most strin- 
gent bounds on variations of fundamental couplings cor- 
related with the gravitational potential [STj- Signifi- 
cant improvements require space-based experiments [38| , 
whether based on clocks or gravitation. 
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Unless otherwise stated, the Einstein summation conven- 
tion will apply to parameter indices i, j, etc.. 
For composite test masses, the differences of the averaged 
values (/p), {Z{Z - 1)^-'*''=') and {A'^'''^) apply. 
We do not consider Cr or Mg since these elements were 
subdominant in test mass composition and are suffi- 
ciently similar to Fe and Si, respectively. 
Very similar results follow from considering only the three 
most sensitive experiments (Be-Ti, Fe-Si02, Al-Pt). 
Similar theoretical issues arise in interpreting atomic 
clock experiments, due to the uncertain dependence of 
nuclear magnetic moments on quark masses. 



